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We examine the photoconductivity of an intrinsic graphene associated with far- and mid-infrared 
irradiation at low temperatures. The model under consideration accounts for the excitation of 
the electron-hole pairs by incident radiation, the interband generation-recombination transitions 
due to thermal radiation, and the intraband energy relaxation due to acoustic phonon scattering. 
The momentum relaxation is assumed to be caused by elastic scattering. The pertinent collision 
integrals are adapted for the case of the massless energy spectrum of carriers that interact with the 
longitudinal acoustic mode and the thermal radiation. It is found that the photoconductivity is 
determined by interplay between weak energy relaxation and generation-recombination processes. 
Due to this the threshold of nonlinear response is fairly low. 

PACS numbers: 73.50.Pz, 73.63.-b, 81.05. Uw 



I. INTRODUCTION 

There are two reasons for unusual transport proper- 
ties of graphene (see Refs. [E [I])' the neutrinolike 
dynamics of carriers, which is described by the Weyl- 
Wallace model, |3] and the specific features of scattering 
processes. In particular, the high efficiency of interband 
optical transitions is associated with a high value of the 
velocity vt,w ~ 10* cm/s characterizing the linear dis- 
persion relations for the graphene valence and conduc- 
tion bands. This is because the matrix element of inter- 
band transitions is proportional to fcvK- On the other 
hand, the coupling of carriers to acoustic phonons ap- 
pears to be weak. As a result, non-equilibrium distri- 
butions of photoexcited carriers can readily be realized 
for the energies smaller than the optical phonon energy 
in the low-temperature region. Due to this, a photore- 
sponse of graphene to far- and mid-infrared (IR) irradi- 
ations should be strong and a low threshold of nonlinear 
response takes place. To the best of our knowledge, the 
effects of non-equilibrium carriers under far- or mid-IR 
excitations is not considered before: both experimental 
and theoretical studies of the transport phenomena in 
graphene are performed at weak dc electric fields (see [1] 
and 0, Q respectively) or under optical excitation (see 
and references therein). 

In this paper, we study the far- and mid-IR photo- 
conductivities of an intrinsic graphene at low tempera- 
tures. As known, [1, 0| the intrinsic graphene exhibits 
a maximum of the dark resistance. Hence, the effect 
of photoconductivity of such a material should be fairly 
strong. Since the concentration of carriers in the intrinsic 
graphene at low temperatures is rather small (less 10^° 
cm~^ at 100 K and lower), one can disregard the intra- 
and interband Coulomb scattering processes. Thus, con- 



sidering the photoresponse of the intrinsic graphene at 
low-temperatures one needs to take into account the fol- 
lowing mechanisms [see Fig. 1(a)]: 

(1) the far- or mid-IR interband photocxcitation; 

(2) the generation-recombination processes due to the in- 
terband transitions caused by thermal radiation; 

(3) the intraband quasi-elastic scattering on acoustic 
phonons; 

(4) the scattering due to a static disorder which is 
an essential mechanism of the momentum relaxation, 
iiii 
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FIG. 1: (a) Schematic view of the interband and intraband 
transitions under consideration, (b) The distribution func- 
tions of thermalized and photoexcited carriers, fl^"""^ and 
A/p (dashed and solid curves, respectively), in the intrinsic 
graphene at low temperatures. Dotted peaks show photoex- 
cited carriers distribution. 

The shape of the non-equilibrium energy distribution 
of carriers is determined by interplay between the radia- 
tive generation-recombination processes and the quasi- 
elastic energy relaxation. Indeed, the corresponding re- 
laxation rates, as we will show below, are of the same 
order at the temperatures about 100 K and lower. Since 
both rates are proportional to the density of states and 
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the radiative transitions are temperature- independent, 
the generation-recombination processes not only deter- 
mine the carrier concentration but also affect the carrier 
energy distribution. Due to an effective recombination of 
low-energy carriers, their energy distribution might be- 
come nonmonotonic exhibiting a peak as shown in Fig. 
1(b). Here we restrict ourselves to the linear (with re- 
spect to pumping intensity) response. The photoconduc- 
tivity appears to be proportional to the concentration of 
the photogenerated carriers for the case of short-range 
momentum relaxation. This concentration as well as the 
photoconductivity decrease with the temperature and the 
energy of the photoexcitation. However, these dependen- 
cies are different if the long-range momentum relaxation 
is essential. 

The paper is organized as follows. The basic equa- 
tions governing the photoconductivity in graphene are 
presented in Sec. H. In Sec. Ill, we consider the car- 
rier distribution for the case of linear (with respect to 
the photoexcitation rate) response. The results of calcu- 
lations of the photoconductivity in this case, including 
numerical estimates, are presented in Sec. IV. The con- 
cluding remarks and discussion of the assumptions used 
are given in Sec. V. In Appendix, the collision integrals 
for the interband relaxation associated with thermal radi- 
ation and for the intraband relaxation caused by acoustic 
phonons are derived. 



II. BASIC EQUATIONS 

We describe the linear response of photoexcited car- 
riers to the dc electric field Eq using the steady state 
kinetic equation for the distribution fmiction Fip: 



dF, 



eEo • ^ = ^ J,{F\lp} + G{F\lp}. 



(1) 



Here I corresponds to conduction {l=+l) or valence (/=- 
1) band, p is the 2D momentum, Jj{F\lp} is the collision 
integral for the jth scattering mechanism, (j = D, LA, 
and R correspond to the static disorder, the acoustic 
phonon scattering, or the radiative-induced interband 
transitions, respectively), and G{F\lp} is the interband 
photogeneration rate. Below we restrict ourselves by the 
weak dc field case, so that Fip ~ Fip + AFip, where Fip 
is the symmetric part of the distribution function under 
the photoexcitation and AF;p oc Eq is the asymmetric 
addition. We evaluate the in-planc isotropic generation 
rate, present the kinetic equation that governs the distri- 
butions functions Fip, and discuss the expression for the 
conductivity considered in Q. 



A. Interband photoexcitation rate 

Within the framework of the Weyl-Wallace model, 
the carrier dynamics under the in-plane ac electric field 



exp{—iujt) + c.c. with the frequency oj is described by 
the zero- field Hamiltonian h^w = vsw{^ • p) with the 
pseudospin Pauli matrix & and the harmonic perturba- 
tion operator: 



(o- • E^)e- 



(2) 



Following the general consideration of interband pho- 
toexcitation (see Sec. 53 in Ref. 8), one can obtain the 
interband photogeneration rate into the Zp-state in the 
form: 



I'p' 



X [5 J {eip - e;/p/ -huj) + 5^ {eip - e//p/ -I- huj)] , (3) 

where the energy conservation law has been introduced 
via the function dj{s) with the phenomenological broad- 
ening energy 7. [§] 

Next, we neglect the in-plane anisotropy of Fip as- 
suming that the dc electric field Eq is sufficiently weak, 
so that one can perform the averaging of the ma- 
trix element in Eq. (3) over the in-plane angle: [lo| 

|(±1, p|iTj,_y| =F liP)l^ = 1/2- As a result, the generation 
rate in Ith band takes form: 



G{F\lp} 



{P-lp - Pip) 



(4) 



with G{F+ip} + G{F_/p} = due to the conservation of 
the net number of the carriers. It is convenient to use in 
the following the electron-hole representation that intro- 
duces the electron (e) and hole {h) distribution functions 
according to the standard replacements [ll| : 



1 - i^-i. 



pt 



fhp- 



(5) 



Considering this substitution, we obtain G{F\ + Ip} 
G{f\p} and G{F\ -l,p} ^ -G{f\p} with the same gen- 
eration rate: 



G{/lri = j( 



evQwE^ 



(1 - /ep - fhp) S.y i2viiwP - Tt-^) , 

(6) 



which is symmetric with respect to the electron-hole re- 
placement, fep ^ fhp- 

Considering the case of relatively weak excitation, we 
assume in the following that fe,hp <C 1 in the vicinity of 
p„ = huj/2v!iw- For such a case, G{f\p) ~ Gp and Gp 
are independent of the carrier distribution: 



Gp = j^bbA 



Sp 



(7) 



Here A{E) = ^5^{E), 6p = '-f/2vaw and vqe is the pho- 
toexcitation frequency. By disregarding the intraband 
photoexcitation of electrons and holes, i.e. the Drude 
mechanism of photoexcitation, we suppose that Sp <ti Puj 
or 7 ^ huj. 
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B. Kinetic equation for fp 

Since the collision integrals [see Eqs. (A. 8) and (A. 12) 
in Appendix A] preserve their form when /ep is replaced 
by ffip, the electron and hole distributions in the intrinsic 
material under the photoexcitation are identical (see Fig. 
1). In addition, the elastic scattering does not affect the 
symmetric distribution due to the energy conservation. 
As a result, the kinetic equation for f^p = fhp = fp takes 
form: 



JiiLAim + Jm{f\p} + Gp = 0. 



(8) 



Here the terms Jrla and Jrh describe the relaxation of 
electrons (holes) caused by the phonon and photon ther- 
mostats with the lattice temperature T. These terms are 
calculated in Appendix. The term associated with the 
interband contribution JsR{f\p} can be presented as 

Jm{f\p} = - Ipf - yl^'-^fl (9) 

where the generation (or recombination) rate, Vp''^^ (or 

Vp'^ '), is expressed via the Planck function N2p/p^^ ~ 

[exp(2p/pi3T) — 1]"^: where p^T = T/v^w is the char- 
acteristic thermal momentum. The rate of spontaneous 

radiative transitions, Vp^^\ is presented [see Eqs. (A. 5) 
and (A. 6) in Appendix] as follows: 

he \ c ) m ~ h ^ ' 

where e is the dielectric permittivity. Here, we have in- 
troduced the characteristic radiative velocity Vr- 

To derive the expression for the in-plane isotropic 
collision integral JblaI/Ip} from (A. 12), it is conve- 
nient to transform the transition probabilities (A. 11) into 
Wpp' = {Wpp, + Wp.p)/2 and AWpp, = Wpp, - Wp,p, 
so that JsLA{f\p} can be presented as: 

JfiLA{/|p} = E^PP'(/p'-/p) (11) 
P' 

-i ^ AT^pp, [(1 - /p)/p, + (1 - /p,)/p] . 
p' 

The energy transfer described by Eq. (A. 11) is small 
because s <C vsw where s is the sound velocity. By using 

also the in-plane averaging \E'(p, p')|p — p'|-^ = (p^+p'^ — 

pp')/2, where the overlap factor ^p', p^ is given by Eq. 

(A. 10), one obtains the transition probabilities: 



AW, 



pp' 



pp 



2v^ 



P«,T^5{P-P') 



(12) 

where D is the deformation potential, Ps is the sheet 
density of graphene and is the normalization area. 



Substituting Eq. (12) into Eq. (11) and integrating by 
parts give us the collision integral in the Fokker-Planck 
form HfH: 



J,LA{f\p} = ^ 



B{qe) 



Mqe) 



( ^ 



p^ dp 

2 



dfp ^ fpjl - fp) 
dp 



\Viiw 



VacP 



PUT 



(13) 



Here, we have introduced the rate of energy relaxation, 
I'p^'^'^^ , and the characteristic velocity Vac [compare to Eq. 
(10)]. 

Since the collision integral (13) is proportional to 
p~^d{. . .}/dp, the integration of the kinetic equation (8) 
over the p-axis yields the following normalization condi- 
tion: 



/ dpppJiiii{f\p} + pp^iyssSp = 0, 
Jo 



(14) 



where pp = 2p/T:h^ViiW is the density of states 
in graphene. Deriving Eq. (14), we have used 
dpppA[{p - Puj)/Sp] ~ Pp^Sp if Sp < pu,. 



C. Conductivity 

The asymmetric parts of electron and hole distribution 
functions have similar form, and A/ep = A/;ip = A/p is 
given by 



A/p = - 



(eE • p) 



.fi(m) 



dfp 
dp 



(15) 



Here 



B(m) 



— VdP'^{plc/K) /Tt- is the momentum re- 
laxation time, where Vd and '^(plc/h) characterize the 
strength of disorder and the quenching of the long-range 
scattering. If a random potential, U^, is characterized 

2 

by the correlation function (C/xt^x') = Ud exp{— [(x — 
x')/Zc]^} with the averaged energy Ud and the correla- 
tion length Ic, ^ one obtains 

-^ = ^S^, *(z) = 9^/1 (z^) , (16) 

Ah Vfiw ^ 

where Ii (x) is the first-order Bessel function of an imag- 
inary argument. 

The conductivity, a, is determined by the standard 
formula [111 



4e viiw f,, 



dfp 
dp 



(17) 



where the nonequilibrium distribution fp is determined 
by Eq. (8). Using the relaxation time Tp"*^ introduced 
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by Eqs. (15) and (16) and integrating by parts in Eq. 
(17), the conductivity can be presented as follows: 



nh Vd 



oo 



energy distribution can be governed by the first order 
equation obtained from Eq. (21). The boundary condi- 
tion assumes the form Afp=p^ — Af^j with 



(18) 



A/. 



(23) 



Here, we have introduced the function ^{z) = 
-5''(2;)/4'(z)2 and take into account ^'(0) = 1/2. In 
the case of short-range scattering, pqtIc/^ ^ 1, for 
the equilibrium conductivity, which corresponds to the 
distribution fp'^'^'^^ = [exp(p/pBT) + 1]"^, one obtains 
cTeq ^ {e'^/Trh)2viiw/vd- 



III. NONEQUILIBRIUM DISTRIBUTION 

Collecting Eqs. (7-9) and (13) together, we arrive at 
the equation for the symmetric distribution function fp: 



By introducing a connection point Pc according to pqt 
Pc Pun we obtain the following slow- varying solution 
for the interval Pc < p < p^'- 



Puj_ 

p 



Pu. , 

— exp 



2pSTVr f dp 



<Jqe 



exp 



2tv / PfiT 

Vqe V P 



Pi 

p 

-P/PIST 



1 g-Pl/PflT 



(24) 



-Plj/pbt 



Mqe 



For the low-energy region, p < pc, we search a solution 
in the form 



dp 



P 



dfp _^ /p(l - fp) 
dp 



PUT 



-{N2p/p^ + l)f;\ +iysEA 



P-Pu, 



Sp 



Afp = exp 



dpi- 



tanh(p/2py) 



0, (19) 



P Put 
{psr/pf 



Sfp 



:^fp^ 



which should be solved with the zero-flow boundary con- 
dition, [dfp/dp + fp{l - fp)/pe,T]p^oo = 0, and the nor- 
malization condition (14). Since the interband photogen- 
eration is centered in the narrow region |p— p^l < 6p, one 
can integrate Eq. (19) over this region. Neglecting an ex- 
ponentially weak flow at p > p^, we consider below the 
uniform equation (19), without the generation term pro- 
portional to vsE- Instead of the generation contribution, 
we use the boundary condition: 



cosh(p/2pfiT)'' ^ 
where dfp is governed by the second order equation 

d'^Sfp Vr/Vq. 



dp^ p2 sinh(p/pBT) 



5^ = 



(25) 



(26) 



with the parameter Vr/vqe ^ 1. Within the WKB ap- 
proximation [l^, we obtain the solution of Eq. (26) in 
the form: 



dfp ^ fpjl - fp) 
dp 



PRT 



P=P^ 



At low pumping, when fp ~ fp^'^^ + Afp and Afp ^ 1 
[Afp GC vse): Eqs. (19) and (20) can be presented as 



v&e5p 



(20) 



Sfp ~ C exp 



dpi 



pi^smh{pi/piiT) 



(27) 



dp [ dp 



dAfp , Afp 



tanh 



PUT 
VrP^ 



p 



Vqe smh{p/pT) 



2pRT 

^fp- 



(21) 



where the normalization constant C is determined from 
the continuity condition at p = Pc which is transformed 
into {piiT/Pc)'^5fp^ = Afp^ cosh(pc/2pflr)- 

As a result, the photoexcited distribution takes the 
form: 



(dAfp Afp 



dp 



P&T J p^p^ 



'^Pi^ Pu) 



Simultaneously, the normalization condition (14) takes 
form: 



'^P \^fp = '^heSpPu. 

cosh[p/ pqt) 



(22) 



In the region p 3> Pbtj where tanh(p/2pj57^) w 1 and 
sinh(p/pi3y) w eP/P<''^/2, one can neglect (dAfp/dp) in 
comparison to Afp/pex- Consequently, a slow tail of the 



A/^ ~ 
where 



exp 



2lir / PBT 

i>,c v P 



-P/P!i 



-Pui/Pe, 



Pc<P < Puj, 

{Aept/[2py cosh(p/2pfiT)]} 
- p 



X exp 



dpi 



exp 



2p 



'i" p Pl\/sinh{pi/psT) 
< p < Pc, 



2«./pbt\ ^.pjp, 

Vqe \ Pc J 



RT 



(29) 



(28) 
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The connection point pc in the distribution (28) is deter- 
mined by 



_ Ac ( Put 


2 V 


Pc 




f 


X exp 










+2 / 




Pc/PBT 




(30) 



Equation (30) is a consequence of Eq. (22) under the 
condition Puj/par ^ 1- The numerical solution of Eq. 
(30) gives the dependency of Pc/pe,T on Vr/vqe which 
can be approximated by the relation Pc/par ~ 0.7 + 
{Vr/Vgef-^/1.5 with an accuracy about 5 %. 



decreases with decreasing of hto and T. For instance, the 
threshold is about 0.6 /xW/cm^ at T =4.2 K and 10 meV 
excitation. 
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FIG. 3: (Color online) Region of parameters (shaded by ho- 
risontal and vertical lines for hu = 100 meV and 40 meV, 
respectively) where the linear regime of response takes place. 
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FIG. 2: (Color online) Distribution function Afp normalized 
to A'' = A/cj(pcj/pfiT)^ under 100 meV photoexcitation at 
T =4.2 K (1), 20 K (2), 40 K (3), and 77 K (4) and under 
10 meV photoexcitation at T =4.2 K (dotted curve). The 
connection points Pc/Pui are marked by arrows. 

The distribution function Afp, which is normalized by 
the factor N = Af^{p^/pQT)'^ and calculated for different 
temperatures and huj = 100 meV, is shown in Fig. 2. 
This function calculated for hoj = 10 meV at T =4.2 K 
is plotted as well. Due to fairly effective recombination 
at p ~ PUT, the distributions approach zero in the low- 
energy region. On the other hand, the fast decreasing 
distributions takes place at p > pc- As a result, the main 
part of photoexcited carriers is localized in the region 
below Pc and a peak of distribution grows up as p^/p^T 
increases. In order to find the conditions of the linear 
(with respect to v^e) response, we use N oc ui^^T^^ 
and set max(A/p) ~ 0.3. Fig. 3 shows the region of 
parameters (in the "pumping - temperature" plane at 
different huj) where the linear response takes place. As 
seen, the threshold of nonlinear response dramatically 



IV. PHOTORESPONSE 

Finally, we consider the photoresponse by using Eq. 
(18) and distribution function (28). The contribution 
of the photoexcited carriers to the conductivity. Act, is 
given by the formula 



.(f)7*PA/, + |/*A/,*(f), (31) 

Pa 

where we have neglected cx Af^^ contribution in the 
lower expression. With the parameter Pujlc/^ = 
huj/{2vBw^/lc), which contains the characteristic energy 
2wfivy^/^c=125 meV at Ic =10 nm, decreases, one can use 
the expansion <i>(z) ~ 4z. In this case Eq. (31) can be 
transformed into the following: 

— ~ 27r;2A„. (32) 

Here An = (4/L^) Afp is the photoexcited concen- 
tration of carriers which is given by 

Pu, 

An ~ -^2 y dppAfp. (33) 



Thus, the photoresponse can be calculated by using the 
substitution of Eq. (28) into Eq. (31) or (33) and nu- 
merical integration. 

Consider first the short-range scattering case, 
Pujlc/Ji *C 1, which is described by the photoinduced 
concentration. An. The temperature dependencies 
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r(K) 

FIG. 4: (Color online) Photoinduced concentration An versus 
temperature: (a) for the pumping intensities 3 /xW/cm^ (1), 
0.3 mW/cm^ (2), and 3 mW/cm^ (3) at ^1^=100 meV and 
(b) for the excitation energies ?ia;=140 meV (1), 100 meV (2), 
and 60 meV (3) at the intensity 1 mW/cm^. Dotted curve 
shows the equilibrium concentration. 



of An, which are calculated for different intensities 
and different frequencies of excitation, are shown in 
Figs. 4(a) and 4(b), respectively. Curves 2 and 3 are 
plotted for the linear response region (see Fig. 3). The 
photoinduced concentration decreases with increasing 
T and fiuj, whereas it increases with pumping intensity. 
Fig. 5 presents the spectral dependencies of An for 
different temperatures. Despite of An exceeds (or 
be comparable with) the equilibrium concentration, 
rieq — 0.52{T/TiVRw)^ 1 the relative photoconductivity 
(32) does not exceed 10^^ for the short-range scattering 
case, when Ic <10 nm for the parameters used. 

In the general case Puilc/^ ^ 1, which occurs if Ic >10 
nm, by using Eqs. (28) and (31) we obtain Aa/aeq vs T 
and huj relations as shown in Figs. 6 and 7, respectively. 
Once again the relative photoconductivity decreases with 
T or tiiu; these dependencies are similar to An vs T and 
tku but an additional dependency on Ic becomes essential. 
In Fig. 7 (curves 2 and 3) we restrict ourselves by the 
linear response region. Since Aa/agq increases with in- 
creasing of Ic or intensity, for the above-discussed region 
of parameters one obtains Aa/acq ^ 1. 

For the case of far-IR excitation with the energy 
huj =10 meV and intensity 1/iW/cm^, one obtains the 
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FIG. 5: (Color online) Spectral dependencies of photoinduced 
concentration for r=50 K (1), 65 K (2), and 77 K (3) at the 
pumping intensity 1 mW/cm^. 
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FIG. 6: (Color online) Relative photoconductivity Aa/a^q vs 
temperature for different correlation lenths lc=W nm (1), 20 
nm (2), and 30 nm (3) at the pumping intensity 3 mW/cm^ 
and huj=100 meV. 

photoinduced concentration An ~ 4.1 • 10^ cm^^ at 4.2 
K while the equilibrium concentration is about 1.7 • 10^ 
cm^'^. Since the short-range scattering condition is valid 
now up to Ic <30 nm, the relative photoconductivity 
Aa/cTcq oc ll does not exceed 0.3 %. 

V. CONCLUDING REMARKS 

In the present work, we have considered the photocon- 
ductivity of intrinsic graphene under far- or mid-IR exci- 
tation of electron-hole pairs. We demonstrated that not 
only carrier concentration but also the energy distribu- 
tion substantially depend on the parameters of excitation 
(frequency and pumping intensity) and the temperature. 
In contrast to the customary semiconductor materials, in 
which the recombination is a most slow process, interplay 
between recombination-generation and energy relaxation 
in graphene is crucial. The carrier distribution smeared 
up to the energy of photoexcitation can be realized due to 
a weak coupling to the phonon thermostat and effective 
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FIG. 7: (Color online) Spectral dependencies of relative pho- 
toconductivity for the pumping intensity 1 mW/cm'^ and the 
correlation lengths lc^2Q nm at r=77 K (1) or r=50 K (2) 
and 20 nm at r=77 K (3) or r=50 K (4). 



(see the discussions in Refs. 4-6) and phonon scatter- 
ing (see [l^ and Refs. therein). Although the analytical 
consideration does not provide a complete qualitative de- 
scription, the relations obtained permit one to estimate 
a character of photoresponsc for the linear regime with 
arbitrary scattering parameters. 

In closing, the obtained results demonstrated a marked 
photosensitivity of intrinsic graphene at low tempera- 
tures. This allows one to analyze the relaxation pro- 
cesses in graphene by using a photoconductivity data. In 
order to check a potential of graphene-based detector ap- 
plications, one needs to perform a numerical modeling, 
including the nonlinear regime of response. 

APPENDIX: COLLISION INTEGRALS 



radiative interband transitions. 

In addition to the physical peculiarities presented, 
the technical results obtained might be useful for the- 
oretical descriptions of the nonequilibrium carriers in 
graphene (heating under dc field, nonlinear optical prop- 
erties, etc). We have evaluated the generation rate that 
describes the interband pumping of electron-hole pairs 
and described by the generation-recombination processes 
caused by the interaction of carriers with the thermal ra- 
diation. We also evaluated the Fokker-Planck collision in- 
tegral governed the quasi-elastic energy relaxation due to 
the deformation interaction between carriers and acous- 
tic phonons. 

Let us briefly discuss the assumptions used in our 
treatment. The main restrictions arise from the low- 
concentration approximation, when one can neglect the 
intra- and interband Coulomb scattering. This approach 
corresponds to an intrinsic graphene at low temperatures 
which is sensitive to photoexcitation due to a high resis- 
tance. We also do not consider a disorder-induced chan- 
nel of generation-recombination and restrict ourselves 
taking into account the simplest scattering mechanisms 
(elastic scattering by static disorder and deformation in- 
teraction with acoustic phonons). These restrictions are 
caused by the lack of data on relaxation process. How- 
ever, the model used allows us to describe important pe- 
culiarities of photoresponsc. Other assumptions are a 
rather standard for the calculations of the transport char- 
acteristics. We used the isotropic energy spectrum both 
for carriers and phonons, and consider the quasi-elastic 
electron-phonon scattering. We neglect an interaction 
with optical phonons because the optical phonon energy 
substantially exceeds T and Tilo. These restrictions are 
valid for low temperatures in the far- and mid-IR spectral 
regions. 

Besides this, we restrict ourselves by the analytical 
consideration of the limiting case presented in Sec. Ill 
and do not perform a numerical solution of Eq. (19) with 
a subsequent integration of (18). This is because the lack 
of a precise data concerns both momentum relaxation 



Below we will evaluate the collision integrals used in 
Sec. II and our consideration is based on the general 
expression Q: 

J{F\a} = [^"'"(1 - Po.)Fo.' - ^"a'(l - FaO^a] , 

(A.l) 

where Fa is the distribution function over a-state. The 
transition probabilities, Waa', are connected by the de- 
tailed balance requirement: Wa>a = exp[— (£„ — ea')/T] 
xWaa' and are determined through 

Q 

X [{Nq + l)S{ea - - hlJq) + Nq5{ea " + tlUJq)] . 

Here Nq is the Planck distribution of the qth boson 
(phonon or photon) mode of frequency cuq and the ma- 
trix clement (a'lxqla) describes the electron-boson inter- 
action. 

1. Radiative-induced transitions 

First, we consider the interaction with thermal radia- 
tion when the secondary-quantized radiation field should 
be substituted into perturbation operator (2) (see Sees. 
20 and 38 in Q). As a result, the operator Xq in (A.2) 
takes the form: 



Xvm^'"iiw\ 77 (o- • e^,BQ) , (A.3) 

y oJsQeV 

where V is the normalization volume, loqq and e,,j3Q are 
the frequency and the polarization vector of the photon 
mode with the 3D wave vector Q and the polarization r;. 
We have also neglected in (A3) the in-plane momentum 
transfer under interband transitions. Calculations of the 
matrix elements (A.3) are performed similar to Sec. HA 
and we obtain 

I (±1,P|X.,BQ |T1,P) P = vlw^^K\ci\' (A.4) 
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while the intraband transitions are forbidden. The av- 
eraging over polarization and direction of photons gives 

After substituting Eq. (A. 4) to general expression 
(A. 2), we obtain the interband probabilities: 



W^+ip,-ip 
W^-ip,+ip 



2{2Trevawf 



3e 



dQ 



NsQ + 1 



(27r)3 '^'y 



(A.5) 



tic phonon scattering while the interband transitions are 
forbidden due to the condition s Wfiiv- The main con- 
tribution to the acoustic phonon scattering appears due 
to the deformation interaction with longitudinal vibra- 
tions [ill, Z5V • Ux, where Ux is the displacement vector 
of LA- mode. By using the quantized displacement oper- 
ator, one obtains the matrix element of Eq. (A. 2) in the 
standard form [11]: 



K^VlXql^P)!' 



5p',p+riq|C'|p-p'|/?i 



(p^p) (A.9) 



where Nqq stands for the Planck number of photons with 
temperature T. Taking the integral over Q-space we ob- 
tain 



W. 



+ip,-ip 



W. 



ip,+ip 



MR) 



2p/pBT 



1 



(A.6) 



Here, q is the 2D wave vector of phonon with frequency 
ujq — sq and |Cqp — D^fiLOq/ {2psS^ L^) stands as the 
electron-phonon matrix element. 

* (p',p), is given by [13: 



The overlap factor, 



where the rate Vp''^^ is given by Eq. (10). 

By using (A.6), one can transform collision integral 
(A.l) as follows 



1 + cos p', p 



(A.IO) 



,B(±) 



p i^p 



(l-^^i,p)F±i 
(A.7) 

where the rates of generation (fp''"*^') and recombination 
i^p^~^) ^'^^ given by Eq. (9). We further use the electron- 
hole representation introduced by Eq. (5) and the radia- 
tive collision integral for electrons takes the form: 

Jm{f\ep} = i^l^~Kl- fep){l~ fhp)~vl^+'^ fephp (A.8) 

and JQR{f\hp} — — JBfl{/|ep}, which are in agreement 
with the particle conservation law. Here the recombina- 
tion term is proportional to fepfhp, while the generation 
contribution is proportional to (1 — /ep)(l — fhp)^ also 



Vp'^^^ and Vp^ ' mean the rates of spontaneous emission 
and absorption. 



Since (A.9) does not depend on I = ±1 (electron-hole 
P' symmetry), we obtain W+ip,+ip' = W^-ip^_ip' = Wpp' 
and the transition probability is given by 

Wpp' = y*(p7p') |qp-p'|A|' (A.ll) 
X {(^IP-P'l/r. + mvQW{p ~ p') - s\p - p'|] 
+N\p-p'\/hS[vswiP - P') + s\p - P'\]} , 

where Nq is the Planck distribution of phonons with tem- 
perature T. By using the electron-hole representation, 
see Eq. (5), one transforms collision integral (A.l) into 
the form: 

JsLA{f\kp} = J2 [^P'P(1 - /'^P)/^?' 



2. Acoustic phonon scattering 



-Wpp,{l-fkp')fkp], (A. 12) 



We further consider the collision integral that is de- 
scribed the intraband transitions caused by the acous- 



with transition probabilities (A.ll) that are the same for 
electrons {k = e) and holes {k = h). 
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